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1.  INTRODUCTION 


The  electrothermal-chemical  (ETC)  propulsion  concept  represents  an  advanced  Aveapon 
technology  in  which  a  conventional  gun  tube  and  combustion  chamber  are  used  to  accelerate  a 
projectile.  The  system  differs,  however,  from  typical  caimon  weaponry  in  that  it  relies  upon  the 
discharge  of  electrical  energy  from  a  power  supply  through  an  insulating  capillary  tube  forming 
an  electrical  plasma,  which  is  injected  into  the  gun  chamber  for  ignition  and  combustion  control  of 
propellant.  The  discharge  forms  an  arc  that  generates  a  hot,  high  pressure,  low  molecular  weight 
plasma  source.  The  high  temperature  plasma  in  the  interior  of  the  capillary  causes  ablation  and 
vaporization  of  surroxmding  insulation,  which  is  enveloped  by  hot  gaseous  material  and  sustains 
the  original  plasma  arc.  The  energetic  plasma  is  then  injected  into  a  bed  of  chemical  propellant  in 
the  gun  combustion  chamber.  The  input  plasma  is  used  to  first  ignite  the  propellant  and  later,  it 
is  proposed  to  drive  and  control  the  combustion  process  of  the  gun.  It  is  hypothesized  that 
control  of  the  combustion  could  result  from  bum  rate  modifications  of  the  propellant  as  a  result 
of  interaction  with  the  high  temperature  plasma.*  The  effects  of  electrical  plasmas  on  bum  rates 
of  some  JA2  sample  propellants  have  been  investigated  in  closed  chamber  ETC  experiments  at 
the  U.S.  Army  Research  Laboratory  (ARE).  The  results,  although  preliminary  and  still  being 
studied,  indicated  potential  bum  rate  increases  compared  to  those  in  conventional  JA2  tests.^ 
Because  of  the  importance  of  the  plasma  source  with  regard  to  the  ETC  technology,  it  was 
decided  to  examine  some  of  the  basic  properties  of  plasma  generation  as  well  as  characteristics  of 
the  plasma  itself,  which  are  described  later  in  this  report. 

1.1  Pulsed  Power  Plasma  Interior  Ballistic  fPPIB)  Code  Description 

PPIB  is  a  time-dependent,  lumped  parameter,  end-to-end  computer  model  of  an  ETC  gun 
system.3  PPIB  simultaneously  models  the  performance  of  the  pulse  power  system,  the  plasma 
cartridge,  and  the  interior  ballistics  of  a  solid  propellant  ETC  gun.  It  is  a  linkage  of  three 
established  models:  (1)  P2SIM,4  (2)  Powell’s  plasma  cartridge  code,^  and  (3)  IBHVG2,6  for 
interior  ballistic  calculations.  The  pulse  power  code  (P2SIM)  is  linked  to  the  plasma  cartridge 
code  by  having  the  plasma  cartridge  act  as  the  pulse  power  system  load.  Any  change  in  the 
output  of  the  pulse  power  system  (e.g.,  increase  in  current)  is  directed  into  the  plasma  cartridge 
and  will  cause  a  change  in  the  state  of  the  plasma.  Any  variation  in  the  plasma  (e.g.,  change  in 
temperature)  will  be  reflected  in  the  pulse  power  system  as  a  change  in  plasma  cartridge 
resistance,  thereby  altering  the  pulse  power  system  load.  Similarly,  any  change  in  the  plasma 
discharge  into  the  gun  will  cause  a  concurrent  change  in  the  combustion  of  the  gun  propellant,  the 
gun  chamber  pressure,  and  ultimately  the  projectile.  A  change  in  the  combustion  of  the  gun 
propellant  is  reflected  back  into  the  plasma  cartridge  as  a  change  in  the  pressure  inside  the 
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cartridge,  which,  during  certain  conditions  (e.g.,  unchoked  flow)  affects  the  output  of  the  cartridge 
and  the  state  of  the  plasma.  This  will  also  affect  the  thermodynamic  properties  of  the  plasma. 
Obviously,  this  change  in  plasma  property  is  reflected  back  to  the  pulse  power  system  as  a 
change  in  the  plasma  cartridge  resistance. 

1.2  Technical  Methods  Employed 

Given  the  plasma  cartridge  length  and  radius,  its  resistance  is  given  by 


in  which  I  is  the  length  of  the  cartridge,  a  is  the  radius,  and  o  is  the  conductivity  of  the  plasma. 
The  conductivity  is  given  by 


in  which  Hq  is  the  density  of  electrons  in  the  plasma,  e  is  the  electron  charge,  rrig  is  the  electron 
mass,  Mqyi  is  the  collision  frequency  for  electrons  and  neutral  atoms,  and  is  the  collision 
fi-equency  for  electrons  and  ions.  The  resistance  value  from  Equation  1  is  the  load  resistance  for 
the  pulse  power  part  of  PPIB.  The  current  density  J  in  the  plasma  is  given  by 


in  which  I  is  the  current  from  the  pulse  power  system  flowing  into  the  plasma  cartridge.  The 
current  density  is  used  in  the  energy  equation  (Equation  3  in  Reference  4)  for  the  plasma 
cartridge.  This  equation  is 


5  ^  2 

—  pco  t/  +  -CO 

dl,  y  2 


'  a  dq  a 


in  which  p,  p,  and  w  are  the  pressure,  density  and  flow  speed;  and  qrs  are  the  longitudinal 
and  radial  heat  fluxes;  U  is  the  internal  energy  of  the  plasma  (per  unit  mass)  and  ^  is  the 
dimensionless  parameter  z/L. 


The  linkage  of  the  P2SIM  and  Powell  codes  was  accomplished  by  using  P2SIM’s  “user- 
defined  load”  capability  with  Powell’s  code  being  the  load.  In  addition,  with  PPIB,  the  user  has 
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the  option  of  including  the  plasma  cartridge  parasitic  inductance,  the  cartridge  lead  resistance,  and 
the  power  cable  resistance. 

To  simulate  the  start  of  the  event,  an  exploding  wire  is  simulated  in  the  model.  The 
expression  for  temperature  of  the  exploding  wire  is  based  on  the  power  in  the  wire  and  the 
assumption  that  the  energy  leaves  the  wire  as  thermal  (black  body)  radiation,  i.e., 

T  =  (fRlkAf^^  (5) 

in  which  ^  is  the  surface  area  of  the  wire  and  k  is  the  Stefan-Boltzmann  constant.  To  simulate 
the  coupling  of  the  radiant  energy  from  the  exploding  wire  to  the  plasma  liner,  the  absorptivity  of 
the  plasma  liner  is  also  included  as  a  model  parameter. 


The  coupling  between  the  Powell  plasma  cartridge  code  and  IBHVG2  is  likewise  very 
straightforward.  The  plasma  energy,  determined  from  the  time  integral  of  input  plasma  power,  is 
added  directly  to  the  energy  balance  equation.  This  equation  is  modified  for  an  ETC  gun  system 
model  and  is  given  as 


X  fij + X  + 

ij  k 


+  J,mkCpk 

k 


+  L 


(6) 


in  which  the  first  sum  is  over  every  surface  j  of  each  propellant  charge  i;  the  second  sum  is  over 
every  species  in  the  plasma.  The  propellant  gas  mass  is  given  by  w^y,  Cy  is  the  specific  heat  at 
constant  volume  for  the  propellant  gases,  and  2ythe  flame  temperature  of  the  propellant.  The 
plasma  species  mass  is  Cp  is  the  specific  heat  of  the  plasma  species,  and  Tp  is  the  plasma 

temperature.  P  is  the  plasma  power  injection  term  and  the  integral  is  from  the  start  of  the  event 
until  the  time,  t.  and  L  are  the  mean  temperature  of  the  system  and  the  energy  loss  terms, 

respectively. 


The  plasma  output  power,  P,  is  comprised  of  three  different  terms.  The  first  is  the  term 
attributable  to  the  thermal  internal  energy  of  the  plasma  and  is  given  by  the  product  of  the 
plasma  specific  internal  energy  at  the  discharge  orifice,  e  (in  J/kg),  and  the  mass  flux  rate,  F  (in 
kg/s),  at  the  discharge  orifice, 

Pj(t)-e(t)F(t)  (7) 

The  second  term  is  the  kinetic  energy  of  the  plasma  and  is  given  as 

P2(t)=  0.5  F(t)(v(t))2  (8) 

in  which  v(t)  is  the  gas  linear  speed. 
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The  third  term  is  the  work  done  against  the  gas  pressure  and  is  given  as 

PS(t)=PoAv(t)  (9) 

with  Pq  being  the  pressure  at  the  discharge  orifice  and^  the  orifice  area.  The  total  power  is  then 

the  sum  of  these  three  terms. 

There  is  no  attempt  to  calculate  the  effects  of  incomplete  mixing  and  transfer  of  energy  to 
the  propellant  combustion  gas.  Therefore,  the  model  will  tend  to  overestimate  the  amount  of 
energy  transfer.  Aside  from  a  complete  treatment  of  the  mixing  of  plasma  and  propellant  gas, 
which  is  beyond  the  scope  of  this  project,  the  only  way  to  accommodate  incomplete  mixing  as 
well  as  to  control  the  contribution  from  the  terms  in  Equations  8  and  9  is  to  introduce  a  set  of 
energy  “coupling”  factors.  There  is  one  factor  for  each  power  term.  A  coupling  factor  of  one  (1) 
denotes  perfect  coupling  of  the  plasma  power  for  that  contribution  to  the  combustion  gas,  while 
a  factor  of  zero  denotes  no  coupling.  The  default  values  in  the  program  are  1,0,  and  0  for 
Equations  7,  8,  and  9,  respectively,  which  indicates  total  transfer  of  internal  plasma  energy  into 
the  gun  chamber  and  no  energy  transferred  from  kinetic  and  work  energy  components. 

The  execution  of  the  code  is  similar  to  that  of  P2SIM.  However,  termination  of  the  event 
will  be  controlled  by  IBHVG2,  which  in  most  cases  occurs  at  projectile  exit.  Data  are  input  to 
the  program  via  three  separate  data  files:  one  for  pulse  power,  one  for  plasma,  and  one  for 
IBHVG2.  The  pulse  power  input  is  virtually  identical  to  the  P2SIM  input  deck,  and  the  plasma 
input  is  virtually  the  same  as  the  original  Powell  code.  The  interior  ballistic  input  deck  is  similar 
to  an  ETC  IBHVG2  deck.  Output  data  from  the  end-to-end  PPIB  code  consist  of  three  data 
files,  one  for  each  subsystem. 

2.  SIMULATIONS 

2.1  Validation  of  P2SIM  (pulsed  power  simulations') 

P2SIM  was  validated  through  comparisons  with  simulations  obtained  from  electronic 
circuit  analysis  software.  For  this  task.  Microcap  III  (Spectrum  Software)  was  used  as  a 
standard  to  predict  the  performance  of  an  electric  gim  type  pulsed  power  supply  discharging  into 
a  fixed  resistive  load.^  Although  the  fixed  resistive  load  does  not  model  the  exact  electrical 
behavior  of  the  electric  gun  plasma  load,  it  does  provide  a  reasonable  approximation  of  the 
electrical  behavior  of  the  power  supply.  The  results  of  simulations  from  P2SIM  and  Microcap 
III  for  a  4-MJ  pulsed  power  supply  configured  from  a  network  consisting  of  capacitors, 
inductors,  switches,  and  other  components  are  given  in  Figure  1,  where  the  two  sets  of  predicted 
curves  are  within  5%  of  perfect  agreement.  Other  validation  exercises  have  been  performed  for 
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other  power  supply  configurations  that  are  not  detailed  here  but  were  performed  under  contract 
to  the  U.S.  Army. 8  Included  among  the  types  of  power  supply  configurations  that  can  be 
modeled  by  PPIB  are  compulsators  and  homopolar  generators  as  well  as  capacitor-based  pulse¬ 
forming  networks  (PFNs)  as  illustrated  in  this  report.  The  input  deck  for  P2SIM  used  to 
generate  the  data  in  Figure  1  is  given  in  Appendix  A,  Part  1. 
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Figure  1.  Validation  of  PPIB  with  Microcap  III  Showing  Delivered  Electrical  Power  and  Energy 
to  a  120-mQ  Resistive  Load. 


2.2  Validation  of  the  Plasma  Model 

To  validate  the  pulse  power  and  plasma  linkage,  experimental  data  were  compared  to 
predictions  from  the  linked  code.  These  data  were  taken  from  test  firings  from  the  ARL  five- 
module  PFN  with  the  capacitors  charged  to  7  kV,  which  was  discharged  into  an  open  air  type 
plasma  capillary.  The  switch  firing  times  were  set  for  a  short  pulse  at  0  s,  40  s,  80  s,  120  s,  and 
180  s.  The  plasma  cartridge  dimensions  used  were  length  10.92  cm  (4.3  inches)  and  inner 
diameter  0.95  cm  (0.375  inch).  The  model  input  included  the  pulse  power  parameters  (capacitor 
voltages  and  switch  times),  the  plasma  cartridge  parameters  (capillary  radius  and  length),  the 
number  of  iterations  per  time  step  and  the  physics  parameters  of  the  plasma  (partition  function 
data,  ionization  potentials,  etc.).  The  PFN  and  plasma  cartridge  input  decks  for  these 
calculations  are  given  in  Appendix  A,  Parts  2  and  3.  For  the  purpose  of  validation,  the  measured 
load  current  was  compared  to  the  current  predicted  by  the  model.  The  agreement  between  the 
two  current  waveforms  is  quite  good,  with  the  model  slightly  over-predicting  the  peak  current 
value  at  an  earlier  time  than  the  actual  data.  The  predicted  peak  current  is  64.3  kA,  occurring  at 
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approximately  237  s,  while  the  measured  peak  is  61.5  kA,  occurring  at  280  s.  The  reason  for  the 
difference  in  peak  values  is  probably  attributable  to  added  resistive  losses  in  the  actual  system, 
and  the  difference  in  time  is  most  likely  attributable  to  inductive  effects  in  either  the  pulsed 
power/plasma  system  or  in  the  measurement  apparatus. 

In  addition,  the  linked  plasma  model  was  further  validated  by  comparing  output  data  with 
the  stand-alone  Powell  plasma  model.  For  these  simulations,  the  model  input  parameters 
adjusted  were  the  current  magnitude,  capillary  radius,  and  capillary  length.  The  predicted  bulk 
plasma  temperature  and  plasma  resistance  from  the  two  models  are  shown  in  Table  1 .  The 
conductivity  model  assumed  for  these  calculations  is  that  of  Kurilenkov  and  Valuev,^  and  a 
surface  temperature  lower  than  the  bulk  plasma  temperature  is  also  assumed.  Note  that  the 
linked  plasma  model  in  PPIB  incorporates  the  entire  end-to-end  code  including  P2SIM  and 
IBHVG2.  Since  it  is  necessary  to  test  the  plasma  model  in  a  condition  independent  from  the 
interior  ballistic  portion  of  the  code,  a  large  chamber  volume  containing  a  small  propellant  charge 
was  modeled  by  the  IBHVG2  subroutine.  This  allowed  for  a  choked  flow  condition  in  the 
plasma  model  for  the  duration  of  the  calculation,  which  was  the  case  for  the  stand-alone  model. 
The  results  of  the  calculations  from  the  two  codes  are  given  in  Table  1,  where  the  current 
magnitude  varies  from  16  kA  to  287  kA.  Once  again,  the  agreement  between  the  two  models  is 
reasonably  good  with  the  largest  discrepancy  of  about  13%  occurring  for  the  resistance  in  the 
150-kA  case.  One  possible  explanation  for  this  discrepancy  is  the  fact  that  the  Powell  model  as 
used  here  provides  output  for  a  single  input  current  data  point,  whereas  the  PPIB  code  predicts 
output  for  the  entire  current  waveform.  The  PPIB  code  is  time  dependent  whereas  the  Powell 
model  used  here  is  a  steady  state  code.  The  output  data  from  PPIB  and  Powell  are  compared  at  a 
single  current  data  point  used  in  the  Powell  code.  This  means  that  PPIB,  in  most  cases,  "will  have 
calculated  plasma  data  for  many  points  before  reaching  the  current  level  assumed  in  the  Powell 
code.  Obviously,  any  time-dependent  changes  in  capillary  conditions  could  account  for  errors  of 
this  size  since  the  output  data  are  generally  a  strong  function  of  the  initial  input  parameters. 

Once  again,  a  sample  input  deck  for  PPIB  is  provided  in  Appendix  A,  Part  4. 

2.3  High  Energy  ETC  Plasma  Calculations 

The  Powell  plasma  model  was  used  to  characterize  the  energy  partitioning  of  the  plasma  as 
a  result  of  a  discharge  from  a  portion  of  a  4-MJ  pulsed  power  supply,  specifically  a  2-MJ 
portion  of  the  full  4-MJ  power  supply.  The  results  from  the  calculations  are  given  in  Table  2 
for  each  of  three  cases  in  which  the  capillary  radius  is  varied  from  1.92  to  7.00  mm  for  a  given 
current  amplitude— in  this  case,  75  kA.  The  plasma  code,  which  is  described  in  great  detail  in 
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Reference  4,  predicts  the  plasma  internal,  kinetic,  and  work  energies  based  on  the  input 
parameters  of  plasma  capillary  radius,  length,  and  input  current  magnitude.  The  code  was  used 
here  to  determine  the  amount  of  energy  contained  in  each  of  the  three  energy  components 
(internal,  kinetic,  work)  for  a  realistic  power  supply  configuration  and  with  realistic  dimensions. 
The  capillary  radius  is  fixed  at  1 1 .84  cm  (sized  for  a  120-mm  tank  cannon  application).  The 
internal  energy  consists  of  the  energy  required  to  vaporize,  dissociate,  and  ionize  the  plasma 
material,  thermal  energy  of  the  plasma,  and  electronic  excitation  energy.  The  kinetic  energy  is 
associated  with  the  movement  of  plasma  through  the  capillary  tube  (1/2  m  v^),  while  the  work 
energy  is  the  energy  against  pressure  in  the  capillary  (P/p).  As  seen  in  the  data  of  Table  2,  for  all 
cases,  the  majority  of  energy  is  contained  in  the  internal  energy  component,  with  the  smallest 
component  being  kinetic  energy.  The  general  trend  is  for  a  decrease  in  internal  energy  and  an 
increase  in  kinetic  and  work  energies  as  the  capillary  radius  is  decreased,  which  seems  reasonable 
since  the  capillary  pressure  rises  rapidly  with  decreasing  capillary  radius.  However,  the  internal 
energy  is  fairly  stable  over  the  range  of  radii  examined  with  only  8.5%  change  from  smallest  to 
largest  capillary. 


Table  1.  Comparison  of  PPIB  and  Powell  Plasma  Model  Simulations  fwith  Kurilenko v- Valuev 
conductivity  model  and  surface  temperature  assumption') 


INPUT  PARAMETERS 

CURRENT  CAPILLARY  CAPILLARY 

(AMPS)  RADIUS  LENGTH 

(mm)  (cm) 

PLASMA  OUTPUT 

POWELL  PPIB 

Temp.  Resistance  Temp.  Resistance 

(K)  (mQ)  (K)  (mQ) 

187,000 

1.92 

11.6 

156,349 

54.0 

164,000 

48 

16,000 

3.17 

7.62 

28,407 

121.0 

31,200 

109 

287,000 

4.75 

11.84 

88,466 

20.2 

96,000 

18 

150,000 

4.75 

11.84 

64,000 

31.2 

69,000 

27 

75,000 

4.75 

11.84 

46,220 

46.0 

51,500 

40 
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Table  2.  Plasma  Energy  Partitioning  for  Three  Test  Cases  With  Varied  Radii 


CASE 

NO. 

CAPILLARY 

RADIUS 

(mm) 

PLASMA 

RESISTANCE* 

(mQ) 

PT.ASMA  ENERGY  PARTITIONING  Inercentl 

INTERNAL  KINETIC  WORK 

ENERGY  ENERGY  ENERGY 

1 

1.92 

138.3 

66.5 

11.2 

22.3 

2 

4.75 

48.6 

73.0 

9.0 

18.0 

3 

7.00 

31.8 

75.0 

8.4 

16.6 

*Plasma  resistance  for  an  input  current  amplitude  of  75  kA  and  capillary  length  of  1 1 .84  cm. 


More  calculations  were  performed  to  approximate  the  electrical  transfer  efficiency  from  the 
PFN  by  calculating  the  electrical  efficiency  for  a  fixed  resistance  load  and  a  modified  PFN 
providing  a  square  current  pulse.  The  power  supply  is  modified  as  a  “transmission  line”  type 
supply  having  inductors  of  a  uniform  value  separating  energy  storage  capacitors. 

This  arrangement  is  used  for  simplicity  and  to  allow  for  a  more  direct  method  of 
determining  its  characteristic  impedance.  The  characteristic  impedance  is  now  proportional  to 
the  square  root  of  the  ratio  of  total  inductance  to  total  capacitance  and  is  calculated  here  at  146 
mQ.  The  resistance  obtained  from  the  plasma  model  (shown  in  Table  2)  is  now  used  by  an 
electronic  circuit  analysis  code  to  predict  the  electrical  transfer  efficiency.  The  resulting  values 
are  given  in  Table  3,  column  3.  This  efficiency  is  now  multiplied  by  the  plasma  partitioning 
percentage  numbers  in  Table  2  to  give  the  last  three  columns  of  numbers  in  Table  3  for  plasma 
energy  delivered  from  the  PFN. 

Finally,  Table  4  provides  a  prediction  from  a  black  body  radiation  model  showing  the 
percentage  of  radiation  belonging  to  each  of  the  types  of  radiation.  ^  ^  Note  that  at  the  high 
internal  plasma  temperatures  predicted,  in  all  cases  the  majority  of  the  energy  is  in  the  ultraviolet 
(UV)  region.  It  is  worth  noting  that  similar  internal  (bulk)  temperatures  (3  eV)  for  capillaries  of 
this  size  and  input  current  level  have  been  predicted  by  other  researchers.  Assuming  that  all 
the  energy  in  the  plasma  is  radiative  except  for  the  kinetic  energy  component,  then  for  Case  1, 
approximately  72.8%  of  the  energy  would  be  in  the  form  of  radiative  energy. 


Table  3.  Percentage  of  delivered  electrical  energy  from  the  pulsed  power  supply  to  plasma 
components  of  internal,  kinetic,  and  work  energies 


CASE 

NO. 

CAPILLARY 

RADIUS 

(mm) 

PERCENT 

ENERGY 

DELIVERED 

FROM  PFN* 

PERCENT  ENERGY  TRANSFERRED  FROM  PFN 

INTERNAL  KINETIC  WORK 

ENERGY  ENERGY  ENERGY 

1 

1.92 

82 

(1,361) 

54.5 

(905) 

9.2 

(153) 

18.3 

(304) 

2 

4.75 

57 

(946) 

41.6 

(691) 

5.1 

(85) 

10.3 

(171) 

7.00 

44 

(730) 

33.0 

(548) 

3.7 

(61) 

13 

(121) 

*For  a  square  current  pulse  from  a  2-MJ  PFN  with  1.66-MJ  stored  electrical  energy  having  a  characteristic 
impedance  of  146  mQ.  The  values  in  parantheses  indicate  the  amount  of  energy,  in  kJ,  obtained  from  the  original 
1.66  MJ  of  stored  energy. 


Table  4.  Black  Body  Radiation  Prediction  for  Three  Plasma  Test  Cases 


CASE 

NO. 

RADIUS 

(mm) 

PLASMA* 

TEMP.  (K,  eV) 

PLASMA  RADIATION  PARTITIONING  fPERCEND 

UF  VISIBLE  INFRARED 

(0  -  0.3  micron)  (0.4  -0.7  micron)  (0.8  -  2  micron) 

1 

1.92 

55,700,  4.8 

97.8  (1,182) 

2.1  (25) 

0.1  (1.2) 

2 

4.75 

34,232,  2.9 

92.2  (795) 

7.2  (62) 

0.6  (5) 

3 

7.00 

27,965,  2.4 

87.5  (585) 

11.4  (76) 

1.1  (7) 

*  Similar  plasma  temperature  results  have  been  obtained  from  other  plasma  researchers  (see  Reference  12).  The 
values  in  parantheses  indicate  the  amount  of  energy,  in  kJ,  obtained  from  the  original  1 .66  MJ  of  stored  energy, 
assuming  only  internal  and  work  energies  contribute  to  radiation. 


While  this  assumption  is  somewhat  arbitrary  and  to  be  taken  only  as  an  approximation,  it 
provides  a  simple  way  of  accounting  for  energy  losses  encountered  in  the  generation  of  the 
plasma.  The  results  from  this  analysis  indicate  that,  based  on  the  approximated  power  supply 
efficiency  and  predicted  internal  plasma  temperature  for  Case  1,  of  the  original  1.66-MJ  store  of 
electrical  energy,  approximately  1.18  MJ  are  converted  to  ultraviolet  energy.  The  remaining  27 
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kJ  delivered  to  the  plasma  are  shared  by  visible  and  infrared  energy  with  about  25  kJ  in  the 
visible. 

In  addition,  experimental  work  for  the  purpose  of  determining  physical  interaction  of 
plasmas  and  propellants  is  progressing  at  ARL.  In  these  studies,  gun  propellant  samples  are 
allowed  to  interact  with  the  radiation  from  high  energy  plasmas,  typical  of  those  studied  in  this 
report.  To  date,  the  experimental  results  indicate  a  rapid  optical  reaction  or  response  of  the 
samples  exposed  to  the  plasma  radiation.  These  tests  also  show  signs  of  increased  propellant 
surface  areas  following  plasma  exposure  and  the  importance  of  propellant  location  and  geometry, 
all  of  which  will  also  require  further  investigations. 

3.  SUMMARY 

An  electrothermal-chemcial  (ETC)  gun  model,  designed  by  Princeton  Combustion  Research 
Laboratories  (PCRL),  for  predicting  the  pulsed  power,  plasma,  and  interior  ballistic  aspects  of  an 
ETC  gun  system  has  been  completed  and  delivered  under  contract  to  ARL.  This  model,  PPIB,  is 
configured  as  an  end-to-end  gun  design  code  incorporating  the  three  main  subsystems  of  an  ETC 
gun,  based  upon  the  well-known  models  of  P2SIM,  the  Powell  plasma  model,  and  IBHVG2. 

This  code  has  been  validated  against  several  sets  of  data,  including  pulsed  power  data  from  a  4- 
MJ  system,  30-mm  ETC  plasma  experiments,  and  the  stand-alone  version  of  the  Powell  plasma 
model. 

The  results  from  the  one-dimensional  Powell  plasma  model  indicate  that  the  capillary 
resistance  is  a  sensitive  function  of  input  current  and  capillary  radius.  For  the  calculations 
performed  here,  the  resistance  varied  by  a  factor  of  4.3  over  the  range  of  radii  used  for  a  fixed 
current  amplitude  of  75  kA.  The  partitioning  of  energy  into  the  various  plasma  constituents  was 
shown  to  be  fairly  insensitive  to  changes  in  radius,  changing  at  most  by  8.5%  over  the  range  of 
radii  used.  However,  for  all  capillary  geometries  and  input  electrical  energies  investigated,  the 
predicted  internal  temperatures  are  quite  high  (between  2.4  and  4.8  eV).  Similar  internal  plasma 
temperatures  for  capillaries  of  this  size  and  input  current  level  have  also  been  predicted  by  other 
plasma  researchers.  ^2  -phe  generation  of  a  high  temperature  plasma  as  seen  in  the  test  cases 
discussed  here  has  the  effect  of  converting  most  of  the  stored  energy  in  the  capacitor  bank  into 
ultraviolet  radiation  before  the  plasma  expansion  process.  The  usefulness  of  the  ultraviolet 
radiation,  how  it  changes  or  re-partitions  during  the  combustion  event,  and  how  it  reacts  with 
materials  such  as  gun  propellants  is  not  clear  at  the  present  time  and  will  require  further  study. 
Also,  experiments  characterizing  the  plasma  propellant  interactions  have  been  performed  at  ARL, 
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which  indicate  rapid  optical  reaction  or  response  of  sample  propellants  exposed  to  plasma 
radiation.  These  tests  show  signs  of  increasing  propellant  surface  area  following  plasma 
exposure  as  well  as  the  importance  of  propellant  location  and  geometry,  all  of  which  will  require 
further  investigations. 
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1.  P2SIM  validation  input  deck 


CONTROL  5uS  2  0 
PULSE  5mSl  10101 

ELEMENT  CAPACITOR  GENE  3  0  1  .99  0  HOO.OuF  OmO  0  0  0  0  0  0  752  446072 

ELEMENT  CAPACITOR  GENE  5  0  1  .99  0  HOO.OuF  OmO  0  0  0  0  0  0  752  446072 

ELEMENT  CAPACITOR  GENE  7  0  1  .99  0  1 400.0uF  OmO  0  0  0  0  0  0  752  446072 

ELEMENT  CAPACITOR  GENE  9  0  1  .99  0  HOO.OuF  OmO  0  0  0  0  0  0  752  446072 

ELEMENT  CAPACITOR  GENE  1 1  0  1  .99  0  HOO.OuF  OmO  0  0  0  0  0  0  752  446072 
ELEMENT  INDUCTOR  GENE  2  3  0  0. 0  120uH  lOmO  0  0  0  47  56000 
ELEMENT  INDUCTOR  GENE  4  5  0  0.  0  80uH  lOmO  0  0  0  47  56000 
ELEMENT  INDUCTOR  GENE  6  7  0  0.  0  60uH  lOmO  0  0  0  47  56000 
ELEMENT  INDUCTOR  GENE  8  9  0  0.  0  43uH  30mO  0  0  0  47  56000 
ELEMENT  INDUCTOR  GENE  10  1 1  0  0.  0  30uH  lOmO  0  0  0  47  56000  ELEMENT 
RESISTOR  GENE  1  0  0  0.  1  I20mO 

ELEMENT  CLS  WITCH  GENE  1  2  OmS  lOmO  0. 1  MO  0  0  0  0  18  44000  ELEMENT 
CLSWITCH  GENE  1  4  0.5mS  lOmO  0.  1  MO  0  0  0  0  1 8  44000 

ELEMENT  CLSWITCH  GENE  1  6  0.75mS  lOmO  0.  IMO  0  0  0  0  1 8  44000 

ELEMENT  CLSWITCH  GENE  1  8  1  .OmS  1  OmO  0.  1  MO  0  0  0  0  1 8  44000 

ELEMENT  CLSWITCH  GENE  1  10  I.2mS  lOmO  0.  IMO  0  0  0  0  18  44000 
ELEMENT  DIODE  GENE  0  3  lOV  45mO  OnH  IMO 

ELEMENT  DIODE  GENE  0  5  lOV  40mO  OnH  IMO  ELEMENT  DIODE  GENE  0  7  lOV 
37mO  OnH  IMO  ELEMENT  DIODE  GENE  0  9  lOV  33mO  OnH  IMO  ELEMENT  DIODE 
GENE  0  1 1  lOV  30mO  OnH  IMO 

ELEMENT  BATTERY  GENE  22KV  lOOOKJ  lOOuO  InH  0  IKA  0  lOKg  1000.  1  ELEMENT 
BATMOD  GENE  0  0  0  1  1  0.  0 
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2.  PPIB  validation  input  deck  for  ARL  five  module  PFN 

CONTROL  5uS  2  0 
PULSE  ImS  1  10  10  1 

ELEMENT  CAPACITOR  GENE  3  0  1  .99  0  1020uF  0.  0.  0.  0.  0.  0.  0. 1.  1000.  1 

ELEMENT  CAPACITOR  GENE  4  0  1  .99  0  488uF  0.  0.  0.  0.  0.  0.  0.  1.  1000.  1 

ELEMENT  CAPACITOR  GENE  6  0  1  .99  0  615uF  0.  0.  0.  0.  0.  0.  0.  1.  1000. 1 

ELEMENT  CAPACITOR  GENE  8  0  1  .99  0  195uF  0.  0.  0.  0.  0.  0.  0.  1.  1000.  1 

ELEMENT  CAPACITOR  GENE  10  0  1  .99  0  400uF  0.  0.  0.  0.  0.  0.  0.  1.  1000.  1  ELEMENT 
INDUCTOR  GENE  1  2  0  0.  0  7uH  0.  0.  0.  0.  10. 1000.  1 

ELEMENT  INDUCTOR  GENE  1  5  0  0.  0  16uH  0.  0.  0.  0.  10.  1000. 1 

ELEMENT  INDUCTOR  GENE  1  7  0  0.  0  28uH  0.  0.  0.  0.  10.  1000. 1 

ELEMENT  INDUCTOR  GENE  1  9  0  0.  0  58uH  0.  0.  0.  0.  10.  1000. 1 

ELEMENT  INDUCTOR  GENE  1  1 1  0  0.  0  15uH  0.  0.  0.  0.  10.  1000.  1  ELEMENT 
RESISTOR  GENE  1  0  0  0.  1  35mO 

ELEMENT  CLSWITCH  GENE  2  3  180uS  lOuO  0.  IMO  0.  0.  0.  0.  10.  1000.  1 
ELEMENT  CLSWITCH  GENE  5  4  120uS  lOuO  0.  IMO  0.  0.  0.  0.  10.  1000.  1 
ELEMENT  CLSWITCH  GENE  7  6  80uS  lOuO  0.  IMO  0.  0.  0.  0.  10.  1000.  1 
ELEMENT  CLSWITCH  GENE  9  8  40uS  lOuO  0.  IMO  0.  0.  0.  0.  10.  1000.  1 
ELEMENT  CLSWITCH  GENE  11  10  OuS  lOuO  0.  IMO  0.  0.  0.  0.  10.  1000.  1 
ELEMENT  DIODE  GENE  0  10  lOV  lOuO  OnH  IMO  0.  0.  0.  0.  1.  1000.  1  ELEMENT 
DIODE  GENE  0  8  lOV  lOuO  OnH  IMO  0.  0.  0.  0. 1. 1000. 1 
ELEMENT  DIODE  GENE  0  6  lOV  lOuO  OnH  IMO  0.  0.  0.  0.  1.  1000.  1 
ELEMENT  DIODE  GENE  0  4  lOV  lOuO  OnH  IMO  0.  0.  0.  0.  1.  1000.  1 
ELEMENT  DIODE  GENE  0  3  lOV  lOuO  OnH  IMO  0.  0.  0.  0.  1.  1000.  1 
ELEMENT  BATTERY  GENE  7KV  lOOOKJ  lOOuO  InH  0  IKA  0  lOKg  1000.  1 
ELEMENT  BATMOD  GENE  0  0  0  1  1  0.  0 
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3.  PPIB  plasma  capillary  input  deck 


4.74E-03,101,10.92E-2,0,1 

400,40 1,1,1  0.5,0. 1 ,0.2, 1 .0E4,273  .,1 .0E-04,.001 ,6. 1 73e7  1 .0, 1  .OOE-5 
1  13  3  2 
2.0  0.00 
1.0  0.00 
3.0  16.4 
5.0  43.5 
5.0  10193.70 
1.0  21648.4 
5.0  33735.2 
9.0  60360.0 
3.0  61982.0 
15.0  64091.0 
3.0  68858.0 
15.0  69715.0 
3.0  70744.0 
9.0  71368.0 
2.0  0.00 
4.0  64.0 
12.0  43024.0 
1.0  0.00 
9.0  52349.0 

0. 

0.9 

1.. 1..1. 

5 

0.  .25  .5  .75  1. 

1.. 0..0. 

C  RI,NZ,XL,IOPT 
C  ITMAX,IPRINT,ICHOKE,ICALC 
C  XI  C,X2C,X1  H,CON,T,XL AM,TCUT,RCUT 
C  PARTITION  FUNCTION  DATA 
CPEX 
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4.  PPIB  plasma  capillary  input  deck  (table  1  data) 

4.75D-03,101,11.84D-2,0,1 

400,40 1,1,1  0.5,0.1 ,0.2, 1 .0D4,273  .,1 .0D-04,.001 ,6. 1 73D7  1 .0, 1  .OD-5 
1  13  3  2 
2.0  0.00 
1.0  0.00 
3.0  16.4 
5.0  43.5 
5.0  10193.70 
1.0  21648.4 
5.0  33735.2 
9.0  60360.0 
3.0  61982.0 
15.0  64091.0 
3.0  68858.0 
15.0  69715.0 
3.0  70744.0 
9.0  71368.0 
2.0  0.00 
4.0  64.0 
12.0  43024.0 
1.0  0.00 
9.0  52349.0 
0. 

0.9 

5 

0.  .25  .5  .75  1. 

l.,l.,l. 

C  RI,NZ,XL,IOPT 
C  ITMAX,IPRINT,ICHOKE,ICALC 
C  XI  C,X2C,X1H,C0N,T,XLAM,TCUT,RCUT 
C  PARTITION  FUNCTION  DATA 
CPEX 
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2  ADMINISTRATOR 

DEFENSE  TECHNICAL  INFO  CENTER 
ATTN  DTICDDA 

8725  JOHN  J  KINGMAN  RD  STE  0944 
FTBELVOIR  VA  22060-6218 

1  COMMANDER 

US  ARMY  MATERIEL  COMMAND 
ATTN  AMCAM 
5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 

1  DIRECTOR 

US  ARMY  RESEARCH  LABORATORY 
ATTN  AMSRLOPSDTA/ 

RECORDS  MANAGEMENT 
2800  POWDER  MILL  RD 
ADELPHIMD  20783-1197 

1  DIRECTOR 

US  ARMY  RESEARCH  LABORATORY 
ATTN  AMSRL  OP  SD  TL/ 

TECHNICAL  LIBRARY 
2800  POWDER  MILL  RD 
ADELPHIMD  207830-1197 

1  DIRECTOR 

US  ARMY  RESEARCH  LABORATORY 
ATTN  AMSRL  OP  SD  TP/ 

TECH  PUBLISHING  BRANCH 
2800  POWDER  MILL  RD 
ADELPHIMD  20783-1197 

1  COMMANDER 
USA  ARDEC 

ATTN  SMCAR  WAH  T  J  DOMEN 

BLDG  62  NORTH 

PCTNY  ARSNL  NJ  07806-5000 

1  COMMANDER 
USA  ARDEC 

ATTN  SMCAR  AEE  BR  K  KLINGAMAN 
BLDG  1501 

PCTNY  ARSNL  NJ  07806-5000 

1  DIRECTOR 
USA  BMD 
ADV  TECH  CTR 
PO  BOX  1500 
HUNTSVILLE  AL  35807 


NO.  OF 

COPIES  ORGANIZATION 

I  CHAIRMAN 

DOD  EXP  SAFETY  BD 
RM  856  C  HOFFMAN  BLDG  1 
2461  EISENHOWER  AVE 
ALEXANDRIA  VA  22331-0600 

1  DA  OFC  OF  PM 

155MM  HWTZR  M109A6  PALADIN 
ATTN  SFAE  AR  HIP  IP 
R  DE  KLEINE 

PCTNY  ARSNL  NJ  07806-5000 

1  COMMANDER 
USA  ARDEC 

ATTN  AMSMC  PBM  E  L  LAIBSON 
PCTNY  ARSNL  NJ  07806-5000 

3  PEO  ARMAMENTS  PM  TMAS 
ATTN  AMCPM  TMA  K  RUSSELL 
AMCPM  TMA  105 
AMCPM  TMA  120 
PCTNY  ARSNL  NJ  07805-5000 

2  DIRECTOR 

BENET  LAB  USA  WATERVLIET  ARNSL 
ATTN  SARWV  RD  G  CARAFANO 
RHASOENBEIN 
WATERVLIET  NY  12189 

2  DIRECTOR 

BENET  LAB  USA  WATERVLIET  ARNSL 
ATTN  SARWV  RD  P  VOTIS 
M  CIPOLLO 

WATERVLIET  NY  12189 

1  COMMANDER 
USA  AMCCOM 

ATTN  AMSMC  IRC  G  COWAN 
SMCAR  ESM  R  W  FORTUNE 
R  ZASTROW 

ROCK  ISLAND  IL  61299-7300 

2  COMMANDER 
USA  AMCCOM 

ATTN  SMCAR  ESM  R  W  FORTUNE 
R  ZASTROW 

ROCK  ISLAND  IL  61299-7300 

1  COMMANDER 
USA  CECOM 
ATTN  ASQNC  ELC  ISL  R 
R&D  TECH  LIB  MYER  CTR 
FT  MONMOUTH  NJ  07703-5301 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  COMMANDANT 
USA  AS 

ATTN  AVIATION  AGENCY 
FT  RUCKER  AL  36360 

1  DIRECTOR 

HQ  TPC\C  RPD 

ATTN  ATCD  MA  MAJ  WILLIAMS 
FT  MONROE  VA  23651-5143 

1  HQ  USAMC 
ATTN  AMCICP  AD 
M  F  FISETTE 

5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 

2  COMMANDER 
USA  ARDEC 

ATTN  SMCAR  CCD  D  SPRING 
SMCAR  CCS 

PCTNY  ARSNL  NJ  07806-5000 

2  COMMANDER 
USA  ARDEC 

ATTN  SMCAR  CCH  T  L  ROSENDORG 
SMCAR  CCH  V  E  FENNELL 
PCTNY  ARSNL  NJ  07806-5000 

2  COMMANDER 
USA  ARDEC 

ATTN  SMCAR  AEE  B  D  DOWNS 
PCHIU 

PCTNY  ARNSL  NJ  07806-5000 

2  COMMANDER 
USA  ARDEC 

ATTN  SMCAR  AEE  J  LANNON 
SMCAR  AED  S  KAPLOWITZ 
PCTNY  ARSNL  NJ  07806-5000 

2  COMMANDER 
USA  ARDEC 

ATTN  SMCAR  FSA  T  M  SALSBURY 
T  GORA 

PCTNY  ARSNL  NJ  07806-5000 

2  COMMANDER 
USA  ARDEC 

ATTN  SMCAR  B  KNUTELSKY 
A  GRAF 

PCTNY  ARSNL  NJ  07806-5000 


1  COMMANDER 
USA  ARDEC 
ATTN  SMCAR 
K  CHEUNG 

PCTNY  ARSNL  NJ  07806-5000 

2  COMMANDER 
USA  ARDEC 

ATTN  SMCAR  EG  G  FERDINAND 
R  ZIMANY 

PCTNY  ARSNL  NJ  07806-5000 

1  US  ARMY  RSCH  LAB 
ATTN  AMSRL  WT  N 
J  INGRAM 

2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1 197 

1  US  ARMY  RSCH  LAB 
ATTNAMSRL  WT  NH 
J  CORRIGAN 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1 145 

1  US  ARMY  RSCH  LAB 
ATTN  AMSRL  WT  NG 
T  OLDHAM 

2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1 145 

1  US  ARMY  RSCH  LAB 
ATTN  AMSRL  EP  MC 
D  SINGH 

BLDG  2707 

FT  MONMOUTH  NJ  07703-5000 

2  COMMANDER 
USARO 

ATTN  TECH  LIB  D  MANN 
PO  BOX  12211 

RSCH  TRI  PK  NC  27709-221 1 

1  COMMANDER 
USABRDEC 

ATTN  STRBE  WC  TECH  LIB  VAULT 
BLDG  315 

FT  BELVOIR  VA  22060-5606 

1  COMMANDER 

USA  TRAC  FT  LEE 
DEFENSE  LOGISTICS  STUDIES 
FT  LEE  VA  23801-6140 
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NO.  OF 

COPIES  OPCiANTZATTON 

1  PRESIDENT 

USA  ARTILLERY  BD 
FT  SILL  OK  73503 

1  COMMANDANT 
USACGSC 

FT  LVNWORTH  KS  66027-5200 

1  COMMANDANT 
USASWS 

ATN  REV  &  TNG  LIT  DIV 
FT  BRAGG  NC  28307 

1  COMMANDER 

RADFORD  ARMY  AMMO  PLANT 
ATTN  SMCRA  QA  HI  LIBRARY 
RADFORD  VA  24141 

1  COMMANDANT 
USAFAS 

ATTN  STSF  TSM  CN 
FT  SILL  OK  73505-5600 

1  DEP  CDR  SDC 

ATN  SFAE  SD  HVL 
D  LIANOS 
PO  BOX  1500 

HUNTSVILLE  AL  35887-8801 

3  COMMANDER 
USA  FSTC 

ATTN  AMSCT  MC  3  S  LEBEAU 
C  BEITER 

220  SEVENTH  ST  NE 
CHARLOTTESVILLE  VA  22901 

1  COMMANDANT 
USA  FACS 
ATTN  ATSF  CO  MW 
B  WILLIS 
FT  SILL  OK  73503 

1  OFC  NAVAL  RSCH 

ATTN  CODE  473  R  S  MILLER 
800  N  QUINCY  ST 
ARLINGTON  VA  22217 

2  COMMANDER 
NSSC 

ATTN  SEA  62R 
SEA  64 

WASHINGTON  DC  20362-5101 


NO.  OF 

COPIES  ORGANIZATION 

1  COMMANDER 
NASC 

ATTN  AIR  954  TECH  LIB 
WASHINGTON  DC  20360 

1  NAVAL  RSCH  LAB 
TECH  LIB 

WASHINGTON  DC  20375 

2  COMMANDER 
HSWC 

ATTN  J  P  CONSAGA 
C  GOTZMER 

SLVR  SPRNG  MD  20902-5000 

2  COMMANDER 
NSWC 

ATTN  CODE  RI3  K  KIM 
CODE  RI3  R  BERNECKER 
SLVR  SPRNG  MD  20902-5000 

3  CDR  NSWC  INDIAN  HEAD  DIV 
ATTN  6210  C  SMITH 

6210J  K  RICE 

62I0C  S  PETERS 

INDIAN  HEAD  MD  20640-5035 

3  CDR  NSWC  DAHLGREN  DIV 
ATTN  CODE  G33  T  DORAN 
J  COPLEY 

CODE  G30  G&M  DIV 
DAHLGREN  VA  2448-5000 

3  CDR  NSWC  DAHLGREN  DIV 
ATTN  CODE  G30I  D  WILSON 
CODE  G32  GUNS  SYSTEMS  DIV 
CODE  E23  TECH  LIB 
DAHLGREN  VA  2448-5000 

1  CDR  NSWC  CRANE  DIV 
ATTN  CODE  4052  S  BACKER 
BLDG  108 

CRANE  IN  47522-5000 

2  CDR  NUSC  ENG  CONV  DEPT 
ATTN  CODE  5B33I  R  S  LAZAR 
TECH  LIB 

NEWPORT  RI  02840 

1  CDR  NSWC  INDIAN  HD  DIV 

ATTN  CODE  270PI  MR  E  CHAN 
101  STRAUS  AVE 
INDIAN  HEAD  MD  20640 
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NO.  OF 

COPIES  ORGANT7ATTON 


NO.  OF 

COPIES  ORGANIZATION 


1  CDR  NSWC  INDIAN  HD  DIV 
ATTN  CODE  3120  MR  R  RAST 
101  STRAUS  AVE 
INDIAN  HEAD  MD  20640 

1  CDR  NSWC  INDIAN  HD  DIV 
ATTN  CODE  210P1  R  SIMMONS 
101  STRAUS  AVE 

INDIAN  HEAD  MD  20640 

2  CDR  NSWC  INDIAN  HD  DIV 
ATTN  CODE  6210  S  BOYLES 
NALMEYDA 

101  STRAUS  AVE 
INDIAN  HEAD  MD  20640 

1  CCDR  NAWC 

ATTN  CODE  3891  A  ATWOOD 
CHINA  LAKE  CA  93555 

1  CDR  USAARDEC 

ATTN  SMCAR  CCH  J  HEDDERICH 
BLDG  1 

PCTNY  ARNSL  NJ  07806-5000 

1  OLAC  PL  TSTL 
ATTN  D  SHIPLETT 
EDWARD  AFB  CA  93523-5000 

10  CIA 

OFC  OF  CENTRAL  REFERENCE 
DISSEMINATION  BRANCH 
RM  GE47  HQS 
WASHINGTON  DC  20502 

1  CIA 

ATTN  J  E  BACKOFEN 
HQ  RM  5F22 
WASHINGTON  DC  20505 

3  DIRECTOR  LANL 
ATTN  B  KASWHIA 
H  DAVIS 

E526  W  REASS 

LOS  ALAMOS  NM  87545 

1  DIRE  LLNL 

ATTN  MS  L355  A  BUCKINGHAM 
PO  BOX  808 
LIVERMORE  CA  94550 


1  DIR  SANDIA  NAT  LAB 

COMBUSTION  RSCH  FACILITY 
ATTN  R  ARMSTRONG 
DIV  8357 

LIVERMORE  CA  94551-0469 

1  DIR  SANDIA  NAT  LAB 

COMBUSTION  RSCH  FACILITY 
ATTNS  VOSEN 
DIV  8357 

LIVEMORE  CA  94551-0469 

1  UNIVERSITY  OF  ILLINOIS 

DEPT  OF  MECH  INDUST  ENGR 
ATTN  PROF  H  KRIER  144  MEB 
1206  N  GREEN  ST 
URB  ANA  IL  61801 

1  JHUCPIA 

ATTN  T  CHRISTIAN 

10630  LTLE  PATUXENT  PKWY 

STE  202 

COLUMBIA  MD  21044-3200 

2  PENN  STATE  UNIV 

DEPT  OF  MECHANICAL  ENGR 
ATTN  J  BROWN 
312  MECHANICAL  ENGR  BLDG 
UNIVERSITY  PK  PA  16802 

1  NCSU 

ATTN  J  G  GILLIGAN 
BOX  7909 

1110  BURLINGTON  ENGR  LABS 
RALEIGH  NC  27695-7909 

1  NCSU 

ATTN  M  BOURHAM 
BOX  7909 

1110  BURLINGTON  ENGR  LABS 
RALEIGH  NC  27695-7909 

2  INST  FOR  ADV  TECH 
ATTN  DR  H  FAIR 

P  SULLIVAN 

4030  2  W  BAKER  LANE 

AUSTIN  TX  78759-5329 

1  INST  FOR  ADV  TECH 
ATTN  DR  I  MCNAB 
4030  2  W  BAKER  LANE 
AUSTIN  TX  78759-5329 
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NO.  OF 

COPIES  ORGANIZATION 

1  SRI  INTERNATIONAL 

PROPULSION  SCIENCE  DIV 
ATTN  TECH  LIB 
333  RAVENSWOOD  AVE 
MENLO  PARK  CA  94025 

1  SPARTA 

ATTN  DR  M  HOLLAND 
945  TOWNE  CENTER  DR 
SAN  DIEGO  CA  92121-1964 

2  UNITED  DEFENSE 
ATTN  MR  M  SEALE 
DRAGIOVANETTI 
4800  E  RIVER  RD 
MINNEAPOLIS  MN  55421-1498 

1  UNITED  DEFENSE 
ATTN  MR  I  DYVIK 
4800  E  RIVER  RD 
MINNEAPOLIS  MN  55421-1498 

1  HERCULES  INC 

RADFORD  ARMY  AMMO  PLANT 
ATTN  D  A  WORRELL 
PO  BOX  1 

RADFORD  VA  24141 

1  HERCULES  INC 

RADFORD  ARMY  AMMO  PLANT 
ATTN  E  SANFORD 
PO  BOX  1 

RADFORD  VA  24141 

2  GDLS 

ATTN  MR  F  LUNSFORD 
DR  M  WEIDNER 
PO  BOX  2074 
WARREN  MI  48090-2074 

2  OLIN  ORDNANCE 

ATTN  V  MCDONALD  LIBRARY 
H  MCELROY 
PO  BOX  222 
ST  MARKS  FL  32355 

1  OLIN  ORDNANCE 

ATTN  D  WORTHINGTON 
PO  BOX  222 
ST  MARKS  FL  32355 

1  PAUL  GOUGH  ASSOC  INC 
ATTN  P  S  GOUGH 
1048  SOUTH  ST 
PORTSMOUTH  NH  03801-5423 


NO.  OF 

COPIES  ORGANIZATION 

1  PHYSICS  INTL  LIBRARY 
ATTN  H  W  WAMPLER 
PO  BOX  5010 

SAN  LEANDRO  CA  94577-0599 

1  ROCKWELL  INTL 
ROCKETDYNE  DIV 
ATTN  BA08  J  E  FLANAGAN 
6633  CANOGA  AVE 
CANOGA  PK  CA  91304 

1  ROCKWELL  INTL 
ROCKETDYNE  DIV 
ATTN  BA08  I  GRAY 
6633  CANOGA  AVE 
CANOGA  PK  CA  91304 

1  SCIENCE  APPLICATIONS  INC 
ATTN  I  BATTEH 
1225  JOHNSON  FERRY  RD 
STE  100 

MARIETTA  GA  30068 

1  SCIENCE  APPLICATIONS  INC 
ATTN  L  THORNHILL 
1225  JOHNSON  FERRY  RD 
STE  100 

MARIETTA  GA  30068 

1  ELI  FREEDMAN  &  ASSOC 
ATTN  E  FREEDMAN 
2411  DIANA  RD 
BALTIMORE  MD  21209 

1  VERITAY  TECH  INC 
ATTN  MR  E  FISHER 
4845  MILLERSPORT  HWY 
E  AMHERST  NY  14051-0305 

1  BATTELLE 

ATTN  TACTEC  LIB 
IN  HUGGINS 
505  KING  AVE 
COLUMBUS  OH  43201-2693 

1  CA  INSTITUTE  OF  TECH 
JET  PROPULSION  LAB 
ATNT  L  D  STRAND  MS  125  224 
4800  OAK  GROVE  DR 
PASADENA  CA  91109 
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1  CA  INSTITUTE  OF  TECH 
JET  PROPULSION  LAB 
ATTN  D  ELLIOT 
4800  OAK  GROVE  DR 
PASADENA  CA  91109 

1  USA  BENET  LAB 

ATTN  SMCAR  CCB  R 
DR  DR  S  SPOCK 
WATERVLIETNY  12189 

1  MARTIN  MARIETTA  DEFENSE  SYS 
ATTN  DR  J  MANDZY 
MAIL  DROP  42  220 
100  PLASTICS  AVE 
PITTSFIELD  MA  01201 

1  STATE  U  OF  NY 

DEPT  OF  ELEC  ENGR 
ATTN  DR  W  J  SARGEANT 
BONNER  HALL  RM  312 
BUFFALO  NY  14260 

1  TEXAS  TECH  UNIVERSITY 

DEPT  OF  EE  COMPUTER  SCIENCE 
ATTN  DR  M  KRISTIANSEN 
LUBBOCK  TX  79409-4439 

1  UNIV  OF  TEXAS  AT  ARLINGTON 
DEPT  OF  ELEC  ENGR 
ATT  DR  L  GORDON 
BOX  19016 

ARLINGTON  TX  76019-0016 

1  WILKES  UNIVERSITY 
DEPT  OF  ELEC  ENGR 
ATTN  DR  A  ARMAND 
WILKES  BARRE  PA  18766 

2  UNIVERSITY  OF  DELAWARE 
DEPT  OF  ELEC  ENGR 
ATTN  DR  N  GALLAGHER 
DR  J  KOLODZEY 
NEWARK  DE  19716 

1  COLORADO  SCHOOL  OF  MINES 
DEPT  OF  ENGINEERING 
ATTN  C  BRAUN 
GOLDEN  CO  80401 

1  SNL 

ATTN  MR  M  GRUBELICH 
DFV  2515 
PO  BOX  5800 

ALBUQUERQUE  NM  87185 


NO.  OF 

COPIES  ORGANIZATION 

2  LOCKHEED  MARTIN  DEF  SYSTEMS 

PRINCETON  COMBUSTION  RSCH  LABS 
ATTN  N  A  MESSINA 
1 1  DEER  PARK  DR  BLDG  IV  SUITE  1 19 
MONMOUTH  JUNCTION  NJ  08852 

1  COMMANDER  USA  ARDEC 
ATTN  SFAEASMTMAAS 

R  BILLINGTON 

PICATINNY  ARSENAL  NJ  07806-5000 

2  RARDE 

GS2  DIVISION  BLDG  R31 

ATTN  DR  C  WOODLEY,  DR  G  COOK 

FORT  HALSTEAD 

SEVENOAKS  KENT  TN14  7BP 

ENGLAND 

1  MATERIALS  RESEARCH  LABORATORY 
SALISBURY  BRANCH 
ATTN  ANNA  WILDEGGER  GAISSMAIER 
EXPLOSIVES  ORDNANCE  DIVISION 
SALISBURY 

SOUTH  AUSTRALIA  5108 
ABERDEEN  PROVING  GROUND 
5  DIRARL 

ATTN  AMSRL  OP  AP  L  (TECH  LIB) 
BLDG  305  APG 

5  CDR  USAATC 
ATTN  WALTON 
GRICE 
D  LACEY 
C  HERUD 
STECS  DA  ID 
L  FRANCIS 

36  DIR  USARL 

ATTN  AMSRL  WT  I  MAY 
J  ROCCHIO 
D  ECCLESHALL 
AMSRL  WT  P  A  HORST 
E  SCHMIDT 
P  KASTE 

AMSRL  WT  PC  BFORCHE 
R  PESCE-RODRIGUEZ 
AMSRL  WT  PB  P  PLOSTINS 
DLYON 
A  ZIELINSKI 

AMSRL  WT  PD  B  BURNS 
AMSRL  WT  PA  T  MINOR 
M  DEL  GUERCIO 
J  DESPIRITO 
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AJUHASZ 
J  KNAPTON 
G  KATULKA  (5  CYS) 

G  KELLER 
P  CONROY 
D  KOOKER 
W  OBERLE 
C  RUTH 

T  ROSENBERGER 
I  STOBIE 
P  TRAN 
K  WHITE 
G  WREN 
M  MCQUAID 

AMSRLWTT  W  MORRISON 
AMSRLWTW  C  MURPHY 
AMSRLWTWB  W  D’AMICO 
AMSRL  WT  WD 
ANIILER 

C  HOLLANDSWORTH 
J  POWELL 


